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We study quadrupole collective oscillations in the bose-fermi mixtures of ultracold atomic 
gases of Yb isotopes, which are realized by Kyoto group. Three kinds of combinations are 
chosen, ^™Yb— ^^^Yb , ^^"Yb— ^^'^Yb and ^^''Yb— ^^"^Yb , where boson-fcrmion interactions 
are weakly repulsive, strongly attractive and strongly repulsive respectively. Collective oscil- 
lations in these mixtures are calculated in a dynamical time-evolution approach formulated 
with the time-dependent Gross-Pitaevskii and the Vlasov equations. The boson oscillations 
are shown to have one collective mode, and the fermions are shown to have the boson-forced 
and two intrinsic modes, which correspond to the inside- and outside-fermion oscillations 
for the boson-distributed regions. In the case of the weak boson-fermion interactions, the 
dynamical calculations are shown to be consistent with the results obtained in the small am- 
plitude approximations as the random phase approximation in early stage of oscillation, but, 
in later stage, these two approaches are shown to give the different results. Also, in the case 
of the strong boson-fermion interactions, discrepancies appear in early stage of oscillation. 
We also analyze these differences in two approaches, and show that they originated in the 
change of the fermion distributions through oscillation. 

PACS numbers: 67.85.Pq,51.10.-Hy 

I. INTRODUCTION 

Over the last several years, there have been significant progresses in ultracold atomic gas physics 
fl: Bose-Einstein condensates (BEG) fl fl 0, 1 

two boson mixtures Q], Fermi-degenerate 



3,Q, 



and so on. In particular, the BE 



atomic gases [3], and Bose-Fermi (BE) mixtures 
mixtures attract physical interests as a typical example where particles obeying different quan- 
tum statistics are intermingled. In the study of this system, we have a big opportunity to 
obtain a lot of new knowledge on quantum many-body systems because we can make a variety 
of mixtures with atomic specie combinations and can control the atomic interactions using the 
Eeshbach-resonance method 11 1|. Theoretical studies of the BE mixtures have been done on 



static properties 
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l],[l3, Q,[3: 



the phase structures and separation 



induced instabilities 



26 



27 



29 



jy the attractive interactions 
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and the collective excitations 



One of important diagnostic signals of many-particle systems is the collective excitations because 
of their sensitivity on the inter-atomic interactions and the ground- and excited-state structures. 
Theoretically, collective motions are usually studied in the random phase ap pro ximation (RPA ) 



261] or its approximate methods: the sum-rule 27|, 123] or the scaling 
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38] 



methods. 

In the previous papers, we calculated time-evolution of the BF mixtures directly with solving the 
time-dependent Gross-Pitaevskii (TDGP) and the Vlasov equations, and studied their monopole 
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321 ] and dipole 33] oscillations; the results are largely consistent with the RPA calculations 



261 ]. but show different behaviors in some aspects: for example, the rapid damping at zero 
temperature (T = 0). 

In the BF mixtures at T = 0, the condensed bosons occupy one single-particle state, and the 



fermions distribute in a wide range of single-particle states. T 
one collective mode in spectrum with no strong damping [39 



lus the boson oscillation has only 



4g]. Especially, 



401 ]: on the other hand, various 

collective modes appear in the fermion collective oscillation generally. The RPA calculations 
{4]] ] can explain the experimental results on the frequencies of collective motions 
in the fermion oscillation of the BF mixtures, two regions of the fermi gas, inside and outside of 
the boson distributions, oscillate with different frequencies, and their interference gives the beat 
and damping phenomena, which are clearly shown in the time-evolution approach; on the other 
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33]. 



hand, the boson oscillation become monotonous 
In the time-evolution approach, the intrinsic frequencies of the oscillation modes can be obtained 
from the Fourier transform of the time-dependence of the collective coordinates. In the dipole 
oscillation of the fermion components, we have really confirmed the existence of the two intrinsic 
modes corresponding to the inside- and outside motions and one forced-oscillation mode caused 
by the boson oscillation [33^. In the early stage of the oscillation, these intrinsic frequencies 
obtained in the time-evolution approach are consistent with the RPA calculation, but, in the 



311 ]: the difference is originated in the 



later stage, the two calculations show different results 
distribution changes of the density and the velocity in time-developments, which appears in 
the time-evolution approach but not in harmonic approximations like the RPA calculations. In 



actual experiments, the oscillation amplitude is not so small 



421 ] that the results obtained in 



the time-dependent approach should be more reliable in comparison with experiments. 
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In this paper, we consider the BF mixtures of the Yb isotopes, which have some particular 
properties; the Yb consists of many kinds of isotopes, five bosons (168, 170,172, 174, lyeyb) g^y^d two 
fermions (^'^^'^''^Yb), which give a variety of combinations in the BF mixtures. Experimental 
researches on the trapped atomic gases of the Yb isotopes are being performed actively by the 
group of Kyoto university; the BEC [43] and the Fermi-degeneracy [4J] have been performed. 



The scatte ring 



the group 
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engths for the boson-fermion interactions have been obtained experimentally by 



45|, 



4a], and the observation of the ground state properties and the collective 
oscillations of the BF mixtures is now under progressing. 

In this paper, we discuss the quadrupole oscillations of the BF mixtures of Yb isotopes in the 
time-evolution approach, where the time evolution of the condensed-boson wave function and 
the fermion phase-space distribution function are obtained from the solutions of the TDGP and 
Vlasov equations, respectively. In the next section, we give the formulation of the transport 
model to calculate the time evolution. In Sec. Ill, the numerical results on the quadrupole oscil- 
lations are shown with their physical properties for three kinds of BF mixtures, ^^''Yb— ^^^Yb , 
i70Yb-i73Yb and I'^'^Yb-i^^yb . Sec. IV is for summary. 



II. TIME-DEPENDENT GROSS-PITAEVSKII AND VLASOV EQUATIONS 



In this section, we briefly explain the time-evolution approach to calculate collective oscillations 
of the BF mixture. Let's consider the system of the coexistent dilute gases of one bosonic and 
one-component fermionic atoms at T = 0, which is trapped in the axially-symmetric potential 
with respect to the z-axis. The zero-range boson-boson and boson-fermion interactions are 
assumed, and no fermion-fermion interaction exists in the system. Then the hamiltonian of the 
system is 



H 



2'KfraBB 
Mb 



2-Kh^aBF{MB + Mf) 



MbMf 



(1) 



where (j) ^-re boson and fermion fields, respectively, Mb,f are the boson and fermion 

masses, il.B,F are the transverse frequencies of the trapping potentials for the boson and the 
fermion, and ubb BF are the boson-boson and boson-fermion s-wave scattering lengths. The 
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transverse and longitudinal components of the spatial coordinate are described by q = (q^., q^), 
and kl is the longitudinal-to-transverse frequency ratio of the trapping potentials; in this paper, 
we assume the same ratio for the bosons and fermions. 

To reduce the parameters, we rewrite Eq. ([T]) with the dimensionless variables; the scaled spatial 
coordinates Rb = {h/MB^sY^'^ and the scaled boson/fermion fields (j) = Rg^^4' & ^ = Rb^^"^^ 
where the scaling parameter Rb is defined by Rb = {h/MB^BY^"^ ■ Then, the scaled hamiltonian 
H = H /HVLb becomes 



H 



+ 



aBB 



{<t>\r)ct>{r)f 



+ hBF<p\r)(j){r)7p\r)ip{r) 



(2) 



where the dimensionless parameters are defined as m/ = Mp/Mb (fermion mass), Uf = VIf/^b 
(fermion trapping-potential frequency), qbb = 8TThaBBR]Y f^BF = ^irhmfaBFi^ + 

mf)~^R'j^ (boson-boson and boson- fermion coupling constants). 

In this paper, we consider the T = system with Ni, bosons and Nj fermions; the total wave 
function of the system is written by 



(3) 



where (pc is the condensed-boson wave function, and is the fermion many-body wave function, 
which is given by the Slater determinant of fermion single-particle wave functions il)n, where n 
is the quantum number. 

The time-evolution equations of the total wave function ^*(t) is obtained from the variational 
condition: 



5 j dT{^{T)\i-^ - H\<^>{t)) = 0. 



(4) 



It gives the coupled TDGP and TDHF equations: 

5 _ . r 1 



{r,T 



1 



2m/ 



(5) 
(6) 
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The effective potentials U b and Up are 

1 

2 
1 



UB{r) = ^{rl + Klz^) + gBBpB{r) + hBFPF{r), (7) 



UfM = ^mfUj){r\ + Klz^) + hBFPB{r), (8) 
where pB and pF are the boson and fermion densities: 

PB{r) = N,\4>c{r)\\ (9) 

occ 

PF{r) = J^lVnWp. (10) 

n 

The number of fermion occupied states in the above equations are usually too large to solve the 
above TDHF equations numerically, so we use the semi-classical approach. In the semi-classical 
limit (/i — > 0), the TDHF equation is proved to be equivalent with the Vlasov equation 

^/(r-,p;r) = ||: + - [V,C/p(r)][Vp]| /(r,p;T) = 0, (11) 

where f{r,p;T) is the fermion phase-space distribution function: 

f{r,p,T)= I d\{^\^{r + ^u,T)^\r-^u,T)\<^)e''P''. (12) 



In actual numerical calculation, we use the test particle method 



48| to solve the Vlasov equation 



pip . In this method, the fermion phase-space distribution function is described by 

f(r,p,T)=^^ J2 S{r-niT)}6{p-p,iT)}, (13) 

^ 1=1 

where Nt is the number of test-particles per fermion. Substituting Eq. (|13p into Eq. (jlip. we 
obtain the equations of motion for the test-particles: 

^n{r) = ^, (14) 

l-p^{r) = -VrUFin). (15) 

Thus the time evolutions of the condensed-boson wave function (f)c and the fermion phase space 
distribution function f{r,p,T) are obtained by solving Eqs. ([5]), (fT4]l and (fT5]l . 



III. GROUND STATES AND TRANSITION STRENGTHS 

Before we show the results on collective excitations in the TDGP+Vlasov approach, we discuss 
the ground and excited states of the BF mixtures of Yb isotopes in RPA calculation. 



6 



The s-wave scattering lengths for the interatomic interactions between Yb isotopes have been 
obtained by Kyoto group {4(]]; these results permit three kinds of BF mixtures, ^^'^Yb— ^^^Yb , 
170yi3_173y)3 and ^''^Yb— ^''^Yb , which are available in the present calculation. They are 
also interesting combinations because the boson-fermion interactions are weakly repulsive 
(i70Yb-i7iYb ), strongly attractive (^^^Yb-^^^Yb ) and strongly repulsive (^^''Yb-^^^Yb ), 
respectively. The values of the coupling constants used in the present calculation are shown in 
TABLE 1. 

For the boson and fermion numbers in the mixtures, we take A*";, = 10000 and A'^^ = 1000, which 
are consistent with the experiments by Kyoto group. In this paper, we assume the spherical 
trapping potentials {hl = 1) with the trapping potential frequency 0,b = 2tt x 100 Hz for bosons 
and fermions w/ = ^Ip/^B = 1- We also use the approximation = Mp/Ms = 1 for all kinds 
of mixtures. 



A. Ground States 

The ground-state wave function of the condensed boson is obtained from the Gross- 
Pitaevskii equation: 

|-^V2 + + gBBPBir) + hBFPF{r)^ #(r-) = Mfe#(r), (16) 

where fii, is the boson chemical potential. To evaluate the fermion phase-space distribution func- 
tion in the ground state, we use the Thomas-Fermi (TF) approximation, where the distribution 
function is assumed to have the form: 

f{r,p) = e[fif-e{r,p)]. (17) 

with the fermion chemical potential /i/. The phase-space one-particle energy e{r,p) in (jl7p is 
defined by 

e{r,p) = -^p^ + UF{r), (18) 
2m/ 

where the effective potential Up is given by ([8]). The fermion density pp included in Up is 
related with the boson density ps through the Thomas-Fermi equation: 

^ {Gn^ppir)}^^^ + ^mjiojirl + kIz^) + hBFPB{r) = Pf. (19) 

Iterating of solving the boson wave function with Eq. ()16p and searching the fermi energy ^/ in 
Eq. (jl9p for the correct fermion number, we obtain the ground-state wave function (p^^^ and the 
fermion density distribution pp to the ground states of the BF mixtures. 
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Thus obtained boson and fermion density distributions are shown in Figs. [T^-c. In the mixtures 
of I'^^Yb-i'^iYb (FigUJi) and ^'^^Yb-^'^^Yb (FigHt), the density distributions are center- 
peaked, but, in the ^^"Yb— ^'^'^Yb mixture, the fermion distribution has large overlap with the 
boson density, because of the attractive boson- fermion interaction (/i_bf < 0). On the other 
hand, in the ^^^Yb— ^^^Yb mixture (FigHb), the fermion density distribution is surface-peaked; 



1 0- 



it is because the BF interaction is strong enough to satisfy the condition hsp/gBB > 
The /iBF-dependence of the fermion density distribution in the TF approximation is explained 
in AppendixEl where the ground-state density profiles are shown to be determined by the single 
parameter hEp/dBB- 

B. Excitation Energies and Transition Strengths in RPA 



As shown in the preceding calculations j27l, l33l|. the collective-oscillation frequency is sensitive 



to the profiles of the density distributions. Here we calculate the excited states in the random 
phase approximation (RPA), which is a good approximation in the case of the small amplitude, 
and show the transition strengths between the excited and ground states. 

For this purpose, we calculate the fermion ground-state wave functions in the Hartree-Fock (HF) 
approximation and, then, the excited states are obtained using the ground state in RPA as in 
Ref. U. 



In the HF approximation, the fermion number Nf is restricted to the values that are determined 
from the subshell closure of fermion single-particle states, so that it become slightly different 
from the value of the present calculation; however, it does not affect the final results. The boson 
and fermion transition amplitudes from an excited state with the excitation energy u;„ to 
the ground state |<I>o) are defined by 

TbK) = ($„|Ob|$o), (20) 
TpiuJn) = ($„|Of|$o>, (21) 

where the quadrupole operators Ob,f are 

Ob = j dM*{r){rl-2z'')(^{r), (22) 

Of = j d^ril)*{r){rl-2z^)^{r). (23) 

(24) 

In Fig. [21 the transition strengths |Tb,fP in RPA are shown for ^'''^Yb— ^''^Yb (left panels), 
170yi3_173y)3 (center panels) and ^^^Yb— ^^'^Yb (right panels) mixtures, where the boson and 



8 



fermion transition strengths are plotted in the top and middle panels, respectively. In the bottom 
panels, the fermion strengths in the single particle processes are added. 

As is clearly shown in the top panels, the boson transition strength \Tb{uj)\^ is almost con- 
centrated on one excited state, whose excitation energy is denoted by lo^q. On the other 
hand, the fermion transition strength |rp(a;)p have two peaks, and distribute between them 
(middle panels). In ^™Yb— ^^^Yb (left panels) and ^'''^Yb— ^''^Yb (center panels), the fermion 
strength |Tp(a')p shows no sharp peaks at the boson excitation energy of Wg; it has almost the 
same structure as the strength in the single particle process (bottom panels). In contrast, the 
^'^^Yb— ^'^^Yb mixture shows a visible strength at lJq (Fig. [2)i), which is very small in the single 
particle process (Fig. O). 

The two peaks of the fermion transition strength in the ^^^Yb— ^^^Yb correspond to the fermion 
intrinsic modes. As shown in Ref. 33], the fermion transition strength of the dipole oscillation 
also shows two peaks, which correspond to the fermion motions in the inside/outside regions 
of the bosons (the inside and outside fermion oscillations). Indeed, the second-peak frequency 
of the quadrupole oscillation is w « 2: the value in the pure harmonic oscillator potential. 
According to Ref. [33], we call these two modes the mode-3 (the inside- fermion oscillation) 



and mode-2 (the outside- fermion oscillation), respectively; then, we define the fermion intrinsic 
frequency uJq as the peak position of the mode-3 in the fermion strength function \Tp{uj)\^ . As 
shown in Fig. [2}i, the ]Tp(ti;)]^ includes an additional mode. From Fig. [2^, we can find that it 
correspond to the oscillation forced by the boson intrinsic oscillation; we call it "the boson-forced 
oscillation" (mode-1). 

It should be noted that, in the ^''''^Yb— ^'''^Yb and ^''"'^Yb— ^'''^Yb mixtures, the fermion strength 
functions have two peaks (except the peak at a; = oJq) but no clear peaks at w ~ 2 corresponding 
to the mode-2 (Figs. [2^,h). It is because one of the two peaks (the stronger peak) are broad in 
these mixtures, and the two peaks are merged into one broad peak. 



IV. QUADRUPOLE OSCILLATIONS OF BOSON-FERMION MIXTURES OF YB 

ISOTOPES 

In this section, we show the results of the time-evolution calculations for the quadrupole oscil- 
lations in the BF mixture. 
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A. Time-Evolution of Boson-Fermion Mixtures and Description of Quadrupole 

Oscillations 

As the initial condition of time-evolution at r = 0, we take the boosted state from the ground- 
state condensed-boson wave function (p^^^ and fermion test particle momentum {Pj-\pz^), which 
are defined by 

^,(r-,r = 0) = exp 1^(4 " 2^^)} #(r), (25) 

and 

Pj^{i) = p^\i) + mfLUfXp rjn, Pz{i) = p[^\i) — 2mfuj fXpz. (26) 

In Eqs. (j25p and (|26p . and Xp are the boost parameters. Then, the boson and fermion 
current densities for the initial state become 

3B{r,T = 0) = ^{VT - 2zz)pf{r), (27) 
j^(r,r = 0) = uOfXF{rT-2zz)pf{r). (28) 

In order to discuss the time-dependence of the quadrupole oscillations, we define the quantities 
Qb,f for boson and fermion: 

B^iB,F)-2Rl{B,F) 

^^'^ = ' ^'^^ 

where 

I^{B,F) = ^/ d^rr^TPBAr), (30) 

Rl{B,F) = ^ ! d^z^psAr), (31) 

and Rq{B) and Rq{F) are the root-mean-square radii of the boson and fermion distributions in 
the ground state. 

In the case of no BF interactions {hsF = 0), the quantities Qb and Qp are proportional to 
ABsin(\/2T) and Ai? sin(2a;/r), respectively; they oscillate monotonously with the periods \/2 
and 2a;/. 

When the amplitude is not small, the oscillations are also not simple and include various modes. 
In order to separate these modes, we use the strength functions Sb,f{'^)-, which is defined as the 
Fourier transform of Qb,f- 

/■*/ 

Sb,f{(^) = / dTQB,FiT)smLL>T (32) 
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where we use the oddness Qb,f{t) = —Qb,f{—t)- It is expected from the initial conditions 
in Eqs. (j25p and (j26p because they lead to Qb,f = and dQB,F/dT 7^ at r = 0. For the 
time-integration in Eq. (j32p . we use the interval < r < 200 unless otherwise noted. 
Finally, we mention how the time-dependence of the collective oscillations obtained in the time- 
evolution method should be compared with the strength function in RPA in comparison. The 
initial boosted state corresponding to the initial condition in Eqs. (j25p and (j26p is given by 

|$(r = 0)) = eH'^sOs+^FO^Ii^^^^ (33) 

where |$o) is the ground state, Ob,f are the quadrupole operators defined in Eqs. ([22|) and 
()23p . and hb,f correspond to the boost parameters. Expanding |$(r = 0)) with respect to kb,f, 
Eq. (133]) becomes 

|$(r = 0)) « \^o) + ^{KBdB + KFdF]\<^o) 

= \^o) + '^J2{KBTB{uJn) + KFTF{0Jn)}\<^n), (34) 
n 

where the transition amplitudes TB^Fi^n) are given in Eqs. (|20p and (j2ip . Thus, using the 

strength functions Tb^f, which is obtained in RPA, the time-dependent functions Qb,f are 
represented by 

Qb{t) = nl/p\ yZ'^Bi^^n) {KBTB{uJn) + KFTF{uJn)} sin{uJnT), (35) 

Qpi^) = AT n2/7^\ {KBTBiuJn) + KF^F (c^n)} sin(a;nT) . (36) 

Also, in RPA, the strength functions become 

^^(a;) = ^TB(a;„) {KBrB(w„) -h KFrp(w„)} F„(t/,ti), (37) 

n 

= ^^^TF{uJn){KBTB{'^n) + KFTF{uJn)}Fn{tf,ti), (38) 
ziV f ^ — ' 

J n 

where 

sin(a; - a;„)t/ - sin(u; - a;„)ti sin (w a;„)t/ - sin (cj -h a;„)ii 
i'n[tf,ti) = —, r — ■ 7 . (39) 

2{u) - UJn) 2[UJ + OJn) 

In the case of extremely small amplitude, the TDGP-|-Vlasov (time-evolution) approach and 
RPA should give the same results for the corresponding boost parameters: kb = A/\/2 and 
kf = Ai?; it has been really confirmed in the case of dipole oscillations [33]. However, in the 
present calculation where the amplitude is not so small, the TDHF and Vlasov approach does 
not give the same results with those in RPA; especially, in the later stage of time-evolution. 
Thus, we use the corresponding values of kb[f) in RPA as to reproduce the oscillations in early 
stage of time-evolution in the TDGP and Vlasov approach. 
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B. Quadrupole Oscillations in i^^Yb i^^Yb 

First, we discuss the quadrupole oscillation in ^^'^Yb— ^''^Yb system, where the boson-fermion 
interaction is weekly repulsive. In Fig. [3l we show the time-dependences of the Qb (the upper 
panel) and Qp (the lower panel) for the initial condition = Xp = 0.1, which corresponds to 
the "in-phase" boson-fermion oscillation. The dashed and solid lines represent the results of the 
TDGP-I-Vlasov and RPA calculations, respectively. Furthermore, we also show the results for 
the initial condition of the "out-of-phase" boson-fermion oscillation: Ab = — Air = 0.1in Fig. [H 
In both cases, the results of the TDGP-|-Vlasov calculation agree well with those in RPA in 
the early stage of time-evolution (r < 20). In RPA, the fermion oscillations show clear beat- 
ing; the amplitudes gradually decrease around r ~ 35 and become larger after that. In the 
TDGP+Vlasov approach, the amplitudes are also damped but does not show clear beat phe- 
nomena. 

Next we calculate the strength functions Sb,f for these oscillations using Eqs. (j32p with inte- 
grating over the range of < r < 200. Figs. [5] and [6] show the Sb,f for the oscillations with 
the in-phase and out-of-phase initial conditions; the boson and fermion strength functions in 
TDGP-I-Vlasov calculation are in the upper left (a) and the lower left (b) panels, and those in 
RPA are in the upper right (c) and the lower right (d) panels. 

In these figures, we find that the boson strength functions Sb have only one sharp peak, which 
is consistent with the monotonous behavior of Qb{t) in Figs. [3] and HI In contrast, the fermion 
strength functions Sp has one negative peak aX uj = and two peaks at a; ~ 1, 87 and ~ 1.96; 
these two peaks are positive in the in-phase initial condition and negative in the out-of-phase 
initial condition. 

These behaviors are very similar to those in the dipole oscillations 33]. The peaks at w ~ lJq 



and at w ~ 1.87 correspond to the boson-forced oscillation (mode-1). and the inside-fermion 
oscillation (mode-3), respectively. In addition, we see a small peak at a; ~ 1.96, which correspond 
to the outside-fermion oscillation (mode-2). In TDGP-I- Vlasov calculations, these peaks have 
broader widths than in RPA, so that the strengths of the outside-fermion oscillation (mode-2) 
in TDGP-I- Vlasov calculation are not so clear (Figs. [3)3 anddJj). 

In order to examine the sign-dependence of the boson-fermion coupling constant /ibf, we 
perform an additional simulation with the attractive boson-fermion coupling constant with 
hBF = —1.9680 (and the same values for the other parameters as in ^'''^Yb— ^''^Yb ). In Figs. [7] 
and [HI we show the time-evolutions Qb,f^'^'^ the strength functions Sb,f with the in-phase 
initial condition. The TDGP-I- Vlasov and RPA calculations give almost the same results. In 
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this case, the fermion oscihation shows no beats in both approaches (Fig. [71 bottom). In the 
strength functions (Fig. [8|), the peak of the fermion oscillation at lo = uJq becomes positive, 
and its heights becomes smaller than that in Fig. [5j We also see a peak at cj ~ 1.99 in Sp in 
the TDGP+Vlasov calculation (Fig. [Hb), which corresponds to the outside- fermion oscillation 
(mode-2), but no clear peaks are found around this frequency in the RPA calculation. In the case 
of attractive boson-fermion interaction, more fermions populate in the boson-distributed region, 
and it makes the outside-fermion oscillation (mode-2) strength small. Thus, the strength of the 
mode-2 becomes very small in the small-amplitude oscillations described in RPA. As the am- 
plitude becomes larger, however, more fermions move to the outside from the boson-distributed 
region, and the strength of the mode-2 becomes larger in TDGP+Vlasov approach. 

C. Quadrupole Oscillations in i^^Yb-i^^yb 

In this subsection, we discuss the quadrupole oscillations in the ^^^Yb— ^^^Yb system, where the 
boson-fermion interaction is strongly attractive. 

In Figs. [9l and [TOl we show the time-dependence of the Qb,f in ^^''Yb— ^^^Yb system for two 
kinds of initial conditions, = Xp = 0.1 (in-phase) and = —Xp = 0.1 (out-of-phase), 
respectively. 

Figs. [9] and [10] clearly show that the TDGP+Vlasov calculation gives very different results 
from RPA in this BF mixture. The boson oscillations in the TDGP+Vlasov calculation have 
slight damping, and shorter periods than those in RPA. Furthermore, Qp decrease in the 
TDGP+Vlasov calculation and become negative when t > 40 (bottom panels). In comparison 
with RPA, the TDGP+Vlasov calculation shows slightly shorter periods in the bose oscillations 
and exhibit slight dampings (top panels). 

In order to examine these results further, we calculate the deformation shifts in the longitudinal 
and transversal directions: Axl = V3Rl/Ro — 1 and Axt = \/3/2Rt/Ro — 1, where Rq^l,t 
are the root-mean-square radii of the ground and the oscillation states, which are defined in 
Eqs. ([30]) and ([3T]) . The time-dependences of Axl,t are plotted in Fig. [11] in the case of the 
in-phase initial condition. In these figures, we can find the increase of the fermion radii in both 
directions; namely the expansion of the fermion gas occurs in the oscillation process. We can 
also find Rl > Rt for the fermion oscillations in this expanding process; from the definition 
of Qp in ([29]) . it is just the origin oi Qp < (Figs. [9] and [lOl bottom). This expansion 
can be understood by the fermion overflow from the boson-distributed region; in the ground 
state, fermions populate largely in the boson-distributed region because of the attractive BF 
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interaction, but, in the excited sate, some of them are pushed into the outside in the course of 
oscillation. Such expansion has already been discussed theoretically in the monopole oscillation 



The strength functions Sb,f corresponding to these oscillations are shown in Figs. 1121 (in-phase) 
and [13] (out-of-phase) . respectively. In RPA calculations (Figs. [T2b.d and[T3b.d). the Sb,f have 
a clear peak at u; ~ 1.4 and a broad peak around uj ~ 2.3, which correspond to the boson-forced 
(mode-1) and the inside-fermion (mode-3) oscillations. No peaks exist, which correspond to the 
outside oscillation (mode-2). It is because, in the ground state, the attractive interaction makes 
a large part of fermions populate in the boson-distributed region so that the fermions that exist 
outside region gives very small contribution in Sb,f- 

The TDGP-I-Vlasov calculation gives similar results for the Sb with those in RPA (Figs.[T2h and 
[T3k): however, the Sp are quite different (Figs. [T2b and[T3b). In the TDGP-|-Vlasov calculation, 
the Sp have no clear peak at w ~ 1.4, which corresponds to the boson-forced oscillation (mode-1), 
but have a broad peak appears around uj ~ 2.3. Also the strength below a; < 1.4 monotonously 
increases as to becomes smaller; these monotonous behaviors in the Sp are also caused by the 
expansion of the fermion gases. 

Furthermore, we see an unclear peak around lo ~ 1.85 in the Sp (Fig. [T2]i). As mentioned 
before, the is assumed to be an odd function for the time parameter r, but the calculated 
strength function does not seem to be the odd function around this peak. In order to clarify the 
situation, we calculate the absolute strength function: 



The results are shown in Fig. [T^.b for the oscillations with the in-phase and the out-of-phase 
initial conditions. Two clear peaks are seen at a; = 1.85 and 2.35 in both cases. To make clear 
the origin of these peaks, we divide the oscillation period into three regions < r < 25 (Fl, 
early stage), 25 < r < 50 (F2, intermediate stage) and 50 < r < 100 (F3, later stage), and 
evaluate the Ap{uj) in these stages. As shown in Fig. \T5\ the peak of the Sp appears around 
UJ = 2.35 in the early stage (panel Fl), but disappear in the latter stage (panel F3). In contrast, 
the peak at lo = 1.85 is seen in all stages. 

Because of the large expansion of the fermion gas, which is shown in Fig. [TTl the quadrupole 
oscillations are supposed to couple with the monopole oscillation mode. For examination, we 
calculate the strength function of the monopole oscillation mode: 



32|. 




(40) 




(41) 
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which are plotted in Fig. [TS]for the three stages: < r < 25 (B1,F1), 25 < r < 50 (B2,F2) and 
50 < r < 100 (B3,F3). The fermion strength function has a peak at a; ~ 2.35 in early stage (Fl) 
but it disappear in later stage (F2,F3); on the other hand, the peak at w ~ 1.85 can be seen in 
all stages. Thus, we can find that the peak at w ~ 1.85 in the Sp of the quadrupole oscillation 
corresponds to the monopole mode, which is excited through the expansion of the fermion gas; 
it explains the asymmetry of the Sp around the peak. 

From the above analysis, we can obtain the following picture on the quadrupole oscillation in the 
170y]-,_173y)3 mixture: 1) The strongly-attractive boson-fermion interaction makes a large part 
of fermions populate in the boson-distributed region in the ground state. 2) In the quadrupole 
oscillation, the oscillating bosons causes the fermion overflow from this region and it caused 
the expansion of the fermion gas. 3) The initial conditions in the present calculation breaks 
the rotation symmetry, and, in the large-amplitude oscillation, the angular momentum becomes 
much larger than 2h. It can excite the monopole mode through the oscillation. As a result, the 
oscillations include both the monopole and quadrupole modes. 4) In the course of the fermion 
gas expansion, the intrinsic modes of the fermion quadrupole oscillation (mode-2 and -3) are 
damped and loses its strength in early stage; instead, the monopole oscillation mode grows up. 
5) In addition the expansion process of the fermion gas is slow and not collective, and then many 
small modes appear below uo < 1.5. 

D. Quadrupole Oscillations in ^^^Yb— ^^■^Yb 

Finally, we discuss the quadrupole oscillations in the ^"^^Yb— ^^'^Yb mixture, where the boson- 
fermion interaction is strongly repulsive and, difi'erent from other two mixtures, the fermion 
density distribution is surface-peaked. 

When hBF > 9bb, the eff'ective potential for fermions has the minimum in the surface of 
the boson-distributed region, which produce the surface-peaked fermion-density distributions 
in the ground state; in the oscillation mode, the fermion gas is considered to oscillate in the 
neighborhood of this minimum. 

Figs. [T7] and [18] show the time-dependence of the Qb,f for the in-phase and out-of-phase initial 
conditions; the dashed and solid lines are for the TDGP+Vlasov and RPA calculations. These 
two calculations are consistent within one cycle of oscillation, but become very different after 
that. In the TDGP+Vlasov calculations, the amplitude of the fermion oscillation shows quick 
damping, and the boson oscillations exhibits beating and small damping in the case of the 
in-phase initial condition (Fig. [T8]l . 
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3l[ and dipole 



The similar damping phenomena have already been appeared in the monopole 
oscillations [33[]; however, in these oscillations, the boson oscillations are stable and simply affects 



the fermion motions through the strong repulsive BF interactions, for quick dampings of the 
fermion oscillations decrease the strength of their intrinsic mode soon. 

In Fig. [T9l we show those oscillations in later stage 60 < r < 120 in the cases of the in-phase (a) 
and the out-of-phase (b) initial conditions. We see that the oscillations keep the almost constant 
amplitudes in this stage, and the bosons and fermions oscillate in out-of-phase with the same 
periods of oscillations. These results shows that the intrinsic fermion mode loses the strength 
in early stage, and only the forced-oscillation mode survives for a long time. 
This mixture prefers the BF out-of-phase oscillation; even if the oscillation starts with the in- 
phase oscillation, the relative phase shift occurs and the oscillation becomes the out-of-phase. 
At the time of the phase change, the oscillation accompanies large diffusion. Thus, the Qp shows 
strong damping, and the Qb also shows small damping and beating throughout the process. 
Next, we show the strength functions Sb,f in Fig. [20] (in-phase initial condition) and in Fig. [21] 
(out-of-phase initial condition). In RPA, the boson strengths concentrate on one sharp peak 
(Fig. [20b and Fig. [2Tb). and the fermion strengths have one sharp peak at the same frequency 
with the Sb and have small peaks with broad widths (Fig. [20l i and Fig. [2Tli): in these results, 
the Sp show no strong intrinsic oscillation modes. 

Here we should note that the boson-boson scattering lengths are very different in ^"^^Yb and 
170 as shown in TABLE I; thus, the boson-boson coupling constant qbb ^^^Yb— ^''^Yb mix- 
ture becomes different from that in the other mixtures. As shown in Appendix [X] in the TF 
approximation, the ground state depends only on the ratio kBp/dBB, and not on qbb directly. 
In order to show this coupling-constant dependence, we calculate the quadrupole oscillation of 
the mixture with qbb = 5.598 x 10^^ and Iibf = 7.420 x 10^^, which is chosen to have the 
same value of the coupling-constant ratio with the ^'''^Yb— ^''^Yb mixture fiBp/gBB = 1-325. 
The strength functions of the quadrupole oscillations are shown in Fig. [22l which should be 
compared with Fig. [20] for the ^"^^Yb— ^'^'^Yb mixture. We find that hight of several peaks are 
different, but their positions are the same in these figures. It confirms that the frequencies of 
the oscillation modes are determined by the ratio kBp/gsB only though peak heights depend 
on two parameters qbb and Ubf- 



16 



E. Fermion Intrinsic Modes and Comparison with Sum-rule Approach 

In order to have some additional discussion about the fermion intrinsic modes, we calculate 
the fermion oscillation when the boson motion is frozen. In Fig. [23l we show the results for 
the mixtures of I'^^Yb-^'^iYb (a), i^OYb-^^^yb (b) and i^^Yb-^^^yb (c); the dampings in 
i70Yb-i73Yb and I'^^Yb-^'^^Yb are found to be much faster than that in ^'^^Yb-^'^^Yb . 
The fermion strength functions Sp of these mixtures are shown in Fig. [251 It is found that the 
Sp show a large positive and a small negative peaks. The positions of the positive peak, which 
are different in these three mixtures, are at the inside-fermion oscillation frequencies (mode-3) 
obtained in the previous calculations. The small peaks are located in 2.0 < uj < 2.1; they are 
supposed to be the outside-fermion oscillations (mode-2) in the previous calculations. 
These peaks have broad widths, and it suggests that the oscillation modes corresponding to them 
have large dampings. Particularly, the outside-fermion mode has very low peak, so that it does 
not play any significant roles in the full calculations in the previous section. As the absolute value 
of the boson-fermion coupling increases, the effective potential for fermions has large deformation 
from the harmonic oscillator shape, and the unharmonic effect from this deformation makes large 
dampings of the fermion intrinsic modes, and the boson-forced oscillation mode (mode-1) have 
larger contribution in the full calculation. We show the /isi?-dependence of the boson and 
fermion intrinsic frequencies in Fig. \25\ which are defined as peak positions of the strength 
functions. 

Now we should give some comments on the results by the sum-rule approach In the sum-rule 
approach, the intrinsic frequencies of the boson and fermion quadrupole oscillations are obtained 

by 



,/4a> + 2.J--^1.«p, (43) 



where 



5 



= TTT^ / dM^-^PpT', (44) 
lOvr^mj J 

Vqp = -^hs, J d rr (45) 

The solid and dashed lines Fig. [25] show the frequencies of the fermion and boson quadrupole 
oscillation obtained in the sum-rule approach. 
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The results of the sum-rule approach well reproduce the TDGP+Vlasov and RPA calculations, 
except the fermion intrinsic frequencies at hsp/gBB = 1-3256. It should be noted that, when 
hsp/gBB > 1, the fermion density becomes surface-peaked, and Vqp < 0, so that, in the sum- 
rule approach, the intrinsic frequencies take minimum values at hsp/gBB ~ 1 as shown in 
Fig. [25l and increase hsp/gBB > 1- 

As mentioned before, when hsp/gBB > 1, the density- and velocity-distribution changes occur 
through the oscillation in the case of the large amplitudes, and the potential minimum for the 
fermion, which exists at the border of the boson-distributed region, also changes its position. 
Thus, it causes the smaller intrinsic frequency of fermion {u!q) than the result obtained in the 
sum-rule approach, because the sum-rule approach assumes the small amplitude oscillations, 
and the density-distribution changes through oscillation are not included. 

Anyway, in the case of the repulsively-strong boson-fermion interaction, the peak of the fermion 
strength function is small and broad, and the contribution from the fermion intrinsic mode is 
not so large except early stage in oscillation. 

V. SUMMARY 

In this paper, we have investigated the collective quadrupole oscillation in three kinds of the BF 
mixtures of the Yb isotopes: ^^^Yb-^^^Yb , I'^^Yb-^'^^yb and ^'^'^Yb-^'^^Yb , where the boson- 
fermion interactions are weakly repulsive, strongly attractive and strongly repulsive. In actual 
numerical calculations, we have obtained the time-evolutions of the oscillating mixtures directly 
using the TDGP and Vlasov equations, and compare the results with the RPA calculations. 
Theoretically, these two approaches predict the same modes of oscillations: the intrinsic (mode- 
2,3) and boson-forced (mode-1) oscillations. Nevertheless, the oscillation behaviors are quite 
different in these two approaches, especially in later stage of oscillation. 

When the boson-fermion interaction is weak, the two approaches give almost the same results in 
early stage, but some difference appears in latter stage. In the mixture with very strong BF in- 
teraction, the difference appears also in the earlier stage; in ^^^Yb— ^^^Yb and ^''^Yb— ^''^Yb mix- 
tures, the two approaches are consistent only in the first one or two periods of the oscillations. 
In the case of the strongly-attractive BF interaction (^'^^Yb— ^'^^Yb ), the fermion overflow from 
the boson-distributed region into the outside causes the fermion gas expansion. On the other 
hand. When the BF interaction is strongly repulsive (^''^Yb— ^''^Yb ), the fermion oscillation 
loses the strength of its intrinsic mode soon, and the fermi gas oscillates with the same period 
of the boson gas in later period. 
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The RPA is available only in the case of the small-amplitude oscillations, because it cannot 
trace the density-distribution changes in the course of time evolution. In actual experiments, 
the amplitude is not so small, and the methods in solving the time-dependence process should 
be proper in comparison with experiments. 

In this paper, we assume the spherical trapping potential with kl = 1, but the actual experiments 
will be done with the largely-deformed potential with kl = ^ /6, for example, in Kyoto group. In 
such cases, the breathing oscillations are not decoupled into the monopole and quadrupole modes 



361 ] but into the longitudinal and transverse oscillation modes. The results of the collective 



oscillations of the BF mixtures in the deformed trapping potential will be discussed in another 
paper [49|]. 



Furthermore, we do not take into account two-body collisions and thermal boson effects [s^. In 
the system Ni, ^ Nj at T = 0, the number of thermal bosons are very small, so that the two- 
body collisions are not expected to play any significant roles in the oscillation processes. In the 
actual experiments, which are performed at very low but T ^ temperatures, thermal bosons 
should give some contributions; the introduction of such effects through two-body collision terms 
into our approach should be done in the future 



511, 



19 



APPENDIX A: GROUND STATE OF BOSON-FERMION MIXTURES IN 
THOMAS-FERMI APPROXIMATION 

In this appendix, we briefly explain the ground state of the BF mixture in the Thomas-Fermi 
(TF) approximation. 

In this approximation, the total energy of the BF mixture is given by 
Et = J 1(4 + 4rl)PB{r) + ^p|(r) 

+ hBFpB{r)pF{r)^, (Al) 

The scaled dimensionless variables and parameters are defined by 

1 hsF m)u)gBB . 
h= 2 ' ^= Q 2 (n,r2,AtLr3), 

,8, ,10^2 o^7, ,8^2 



= PB. ep = fiF. (A2) 



The effective boson and fermion numbers are defined by 



j^20^25 6 

Nb = I d'xuB = J. i Nb, (A3) 



35,10 



Finally, the scaled total energy becomes 



2KLmfufg%s ^ 
= Et 

= J d^x i^x^UB + -^n% + ^n]p + x'^np + huBTiF^ , (A5) 

where = \x\^. 

Let's introduce the particle-number constraints into the total energy as E' = ET—eBNB—epNF, 
where the Lagrange multipliers cb and ep are the scaled boson and fermion chemical potentials. 
The variations of E', SE'/Sub = and SE'/Snp = 0, gives the TF equations for the ground-state 
densities nB,p 

UB + hnp = cb — x^, (A6) 

2 

Up + hnB = ep — x^, (A7) 
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The second order variations give the stabihty condition of the TF ground states: 

2 



(5n^ Sn'jp \6nBSnF I 



> 0, (A8) 



which leads to 



< ^. (A9) 

It should be noted that the scaled equations include three parameters, {eB,eF,h) though four 
parameters {fJ-B, fJ-F, 933, hsp) exist in the original unsealed system. 
Using sp = [nF{x)]^^^ and eliminating ub in Eqs. I\A6\ IA7|) . we obtain 

fisp) = 4 - ^^4 - (eF - hcB) + (1 - h)x'^ = 0. (AlO) 
From the derivative function of (jAlOh : 

^=2.,-3.V, = 3.V,(|,-,,), (All) 

we find that df /dsp > for sp < 2 /?>h? ; it is exactly equivalent to the stability condition in ()A9P . 
It means the existence of the solution with the positive value in ()A10p when /(O) < < /(2/3/i^). 
Now we consider the condition that the fermion density has a maximum peak at the surface. 
Differentiate Eq. (jAlOP with respect to x, we obtain 

^^'"^ (3^ " ^ = "^^^ " ^^^'^^ 
Using the stability condition in ()A9p . Eq. ()A12p gives dsp/dx < when h < 1 (case 1), and 
dsp/dx > when h > 1 (case 2); in case 1, the TF fermion density np should have the maximum 
at X = (the center-peaked profile), and, in case 2, x = should be the minimum of the np, 
which should have a maximum outside the boson occupation region, (the surface-peaked profile). 
An extreme case of the surface-peaked fermion densities (h > 1) is the shell-structure profile, 
where fermions are pushed outside and no fermions exist in the central region. The shell- 
structure profile should appear when the TF solution satisfies np{0) < 0; using Eq. (|A10p . we 
find that it occurs when hcB > ep > cb- It should be noted that, when es > ep, there are no 
solutions in Eq. (jAlOp . 

When h < 1 (the center-peaked fermion profile), the fermion density at the center is required to 
be positive, np{0) > 0; otherwise there are no solutions in Eq. (jAlOp . The condition, np{0) > 0, 
gives the restriction of the parameters as follows: 

4 

heB < ep < hcB + i^^- (A13) 
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When < h < 1, there exists the case that the boson density also shows the surface-peaked 
profile. In order to obtain the second derivative of at x = 0, we two-time differentiate 
Eqs. ()A6l IA7p with respect to x: 



d'^UB 



2, , 1 d'^np , 
5[nH0)l-^l 



+ h ^ r, \x=o 



x=0 

:=0 + h 



dx"^ 



-2, 
-2. 



(A14) 
(A15) 



Solving the above equations for d'^UB / dx"^ , we obtain 

^U=o = {-2 + 3/i[n^(0)]i}|l-^/i2[^^(0)]i| \ (A16) 

It shows that, when 2/3/i < sf(0) = [nir(O)]^/^ < 2/2,h?, d'^UB/dx'^ > at x = 0; the boson 
density have the surface-peaked profile (It should be noted that h satisfies 2/3h < 2/3h? because 
we consider the case of < h < 1). The condition that the solution of Eq. (jAlOp satisfies 
2/3/1 < sf{0) < 2/3/1^ is obtained by f{2/3h) < < f{2/3h^) at x = 0; using Eq. (IATO]) . we 
obtain the condition that the boson density has the surface-peaked profile: 

4 



(A17) 



Also, 71^(0) < gives the boson shell-structure condition: cb < hep . 



[1] S. Giorgini, L.P. Pitaevskii and S. Stringari, Rev. Mod. Phys. 80, 1215 (2008) and reference therein. 
[2] E.A. Cornell and C.E. Wieman, Rev. Mod. Phys. 74, 875 (2002); 

W. Ketterle, Rev. Mod. Phys. 74, 1131 (2002). 
[3] F. Dalfovo, et al., Rev. Mod. Phys. 71, 463 (1999). 

[4] C.J. Pethik and H. Smith, "Bose-Einstein Condensation in Dilute Gases", Cambridge University 
Press (2002). 

[5] J.O. Andersen, Rev. Mod. Phys. 76, 599 (2004). 

[6] A.G. Truscott, K.E. Streker, W.l. McAlexander, G.B. Partridge and R.G. Hulet, Science 291, 2570 
(2001); 

G. Roati, F. RiboH, G. Modugno and M. Inguscio, Phys. Rev. Lett. 89, 150403 (2002). 
[7] B. DeMarco and D.S. Jin, Science 285, 1703 (1999); 

S.R. Granade, M.E. Gehm, K.M. OHara, J.E. Thomas, Phys. Rev. Lett. 88, 120405 (2002). 
[8] F. Schreck, et al., Phys. Rev. Lett. 87, 080403 (2001). 

[9] A.G. Tuscott, K.E. Strecker, W.I. McAlexander, G.B. Parridge, and R.G. HuUet, Science 291, 2570 
(2001); 

Z. Hadzibabic, et al., Phys. Rev. Lett. 88, 160401 (2002); ibid 91, 160401 (2003). 



22 



[10] M. Modugno, F. Ferlaino, F. Riboli, G. Roati, G. Modugno, M. Inguscio, Science 297, 2240 (2002); 

Phys. Rev. A 68, 043626 (2003). 
11] H. Feshbach, Ann. Phys. (NY) 19, 287 (1962). 
12] K. M0lnier, Phys. Rev. Lett. 80, 1804 (1998). 

13] M. Amoruso, A. Minguzzi, S. Stringari, M. P. Tosi and L. Vichi, Eur. Phys. J. D 4, 261 (1998). 

14] T. Miyakawa, K. Oda, T. Suzuki and H. Yabu, J. Phys. Soc. Japan 69, 2779 (2000). 

15] M.J. Bijlsma, B.A. Heringa and H.T.C. Stoof, Phys. Rev. A 61, 053601 (2000). 

16] L. Vichi, M. Inguscio, S. Stringari, and G.M Tino, Eur. Phys. J. D 11, 335 (2000). 

17] L. Vichi, M. Inguscio, S. Stringari and G.M. Tino, J. Phys. B: At. Mol. Opt. Phys. 31 L899 (1998). 

18] N. Nygaard and K. M0lnier, Phys. Rev. A 59, 2974 (1999). 

19] X.X. Yi and CP. Sun, Phys. Rev. A 64, 043608 (2001). 

20] L. Viverit, C.J. Pethick and H. Smith, Phys. Rev. A 61, 053605 (2000). 

21] P. Capuzzi and E.S. Hernandez, Phys. Rev. A 66, 035602 (2002). 

22] P. Capuzzi, A. Minguzzi and M.P. Tosi, Phys. Rev. A 68, 033605 (2003). 

23] T. Miyakawa, T. Suzuki and H. Yabu, Phys. Rev. A 64, 033611 (2001). 

24] R. Roth and H. Feldmeier, Phys. Rev. A 65, 021603(R) (2002). 

25] P. Capuzzi and E. S. Hernandez, Phys. Rev. A 64, 043607 (2001). 

26] T. Sogo, T. Miyakawa, T. Suzuki and H. Yabu, Phys. Rev. A66, 013618 (2002); 

T. Sogo, T. Suzuki and H. Yabu, Phys. Rev. A68, 063607 (2003). 
27] T. Miyakawa, K. Oda, T. Suzuki and H. Yabu, J. Phys. Soc. Jpn 69, 2779 (2000). 
28] T. Miyakawa, T. Suzuki and H. Yabu, Phys. Rev. A 62, 063613 (2000). 
29] A. Minguzzi and M.P. Tosi, Phys. Lett. A 268, 142 (2000). 
30] S.K. Yip, Phys. Rev. A64, 023609 (2001). 

31] T. Maruyama, H. Yabu and T. Suzuki, Phys. Rev. A72, 013609 (2005). 

32] T. Maruyama, H. Yabu and T. Suzuki, Lazer Physics, 15, 656 (2005). 

33] T. Maruyama and G.F. Bertsch, Phys. Rev. A77, 063611 (2008) 

34] X.J. Liu and H. Hu, Phys. Rev. A67, 023613 (2003). 

35] T. Maruyama and G.F. Bertsch, Phys. Rev. A73, 013610 (2006). 

36] T. Maruyama and T. Nishimura, Phys. Rev. A75, 033611 (2007). 

37] G.F. Bertsch, Nucl. Phys. A249, 253 (1975); 

D.M. Brink and Leobardi, Nucl. Phys. A258, 285 (1976). 
38] G.F. Bertsch and K. Strieker, Phys. Rev. CIS, 1312 (1976); 

T. Suzuki, Prog. Theor. Phys. 64, 1627 (1980). 
39] F. Chevy, V. Bretin, P. Rosenbusch, K.W. Madison and J. DaHbard, Phys. Rev. Lett. 88, 250402 
(2002). 

40] M. Tohyama, Phys. Rev. A71, 043613 (2005). 
41] T. Kimura, H. Saito and M. Ueda, J. Phys. Soc. Jpn. 68 (1999); 
R. Graham and D. Walls, Phys. Rev. A57, 484 (1998); 



23 



D. Gorddon and CM. Sarvage, Phys. Rev 85, 1440 (1999). 
[42] D.S Jin, J.R. Ensher, M.R. Matthews, C.E. Wieman and E.A. Cornell, Phys. Rev. Lett. 77, 420 
(1996); 

M.-O. Mewes et al., Phys. Rev. Lett. 77, 988 (1996). 
[43] Y. Takasu, et al., Phys. Rev. Lett. 91, 040404 (2003). 

[44] T. Fukuhara, S. Sugawa, and Y. Takahashi, Phys. Rev. A76, 051604(R) (2007). 

[45] T. Fukuhara, Y. Takasu, M. Kumakura and Y. Takahashi, Phys. Rev. Lett. 98, 030401 (2007). 

[46] K. Enomoto, M. Kitagawa, K. Kasa, S. Tojo and Y. Takahashi, Phys. Rev. Lett. 98, 203201 (2007); 

M. Kitagawa, et al., Phys. Rev. A77, 012719 (2008). 
[47] L.P. Kadanoff and G. Baym, "Quantum Statistical Mechanics" (1962), New York. 
[48] R.W. Hockney AND J.W. Eastwood, "Computer simulations using particles", McGraw-Hill, New 

York, 1981; 

C.Y. Wong, Phys. Rev. C25, 1460 (1982). 
[49] T. Maruyama and H. Yabu, in preparation. 
[50] B. Jackson and E. Zaremba, Phys. Rev. A66, 033606 (2002). 
[51] G.F. Bertsch and S. Das Gupta, Phys. Rep. 160 (1988) 189. 
[52] T. Maruyama, W. Gassing, U. Mosel, K. Weber, Nucl. Phys. A573, 653 (1994). 



24 





System 


aBB (nm) (lOOgss) 


aBF (nm) {IQQhBp) 


hBr/gsB 


(1) 


iTOYb - "lYb 


3.4353 ( 5.5976 ) 


1.9680 ( 3.2067) 


0.573 


(2) 


170Yb - "3Yb 


3.4353 ( 5.5976 ) 


-4.3730 (-7.1255) 


-1.273 


(3) 


174Yb _ 173Yb 


5.4630 ( 8.9016 ) 


7.2410 (11.7994) 


1.325 



TABLE I: Scattering lengths (coupling constant) of the BF mixtures of the Yb isotopes. aBB.BF 
{gBB^hsp) are the boson-boson and boson-fermion scattering lengths (coupling constants). 




FIG. 1: The ground-state density distributions of the BF mixtures: (a) i^^Yb-i^^Yb, (b) ^^"Yb-i^^Yb 
and (c) ^''''Yb— ^''•^Yb (c). The dashed and sohd lines represent the density of the bosons and fermions, 
respectively. 
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FIG. 2: The transition strengths \Ab,f\^ in RPA: i^°Yb-i^iYb (left panels), i^^Yb-^^^yb (center 
panels) and ^^'*Yb— ^^'^Yb (right panels). The top, middle and bottom panels are for the boson, fermion 
and the single-particle fermion strengths 
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FIG. 3: Time evolution of the quadrupole oscillations in the ^™Yb— ^^^Yb mixture with the in-phase 
initial condition A_b=Af = 0.1. The upper and lower panels are Qb.f for the bosons and fermions, and 
the solid and dashed lines represent the results of TDGP + Vlasov and RPA calculations. 




FIG. 4: Time evolution of the quadrupole oscillations in the ^™Yb— ^^^Yb mixture with the out-of-phase 
initial condition As — —Xf — 0.1. The panels and the lines are the same as in Fig. [3l 
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FIG. 5: The strength functions of the quadrupole oscillations with the in-phase initial condition in 
i70Yi3_i7iYb mixture (corresponding to Fig. [3]); the boson and ferniion oscillations in TDGP+Vlasov 
calculation (a and b), and in RPA calculation (c and d). 



20 


_ 1 1 1 


1 1 1 1 1 1 1 1 _ 




; (a) 


TDGP + Vlasov : 


15 






10 






5 











1 ' 



(c) 








1 1 

f| pjii-'' 


- (d) 

-.1,1 


1 , 


, 1 , - 



1.0 1.2 1.4 1.6 1.8 2.0 2.2 2.4.0 1.2 1.4 1.6 1.8 2.0 2.2 2.4 

m/Q.g CO / Q.B 

FIG. 6: The strength functions of the quadrupole oscillations with the out-of-phase initial condition in 
i70Yb_i7iYb mixture (corresponding to Fig. 2]); The panels and the lines are the same with Fig. [51 
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FIG. 7: Time evolution of the quadrupole oscillations with the in-phase initial condition in the BF 
mixture with hsF = —1.9680; other parameters are the same as in ^™Yb— ^^^Yb case. The panels and 
the lines are the same as in Fig. \E\ 
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FIG. 8: The strength functions of the quadrupole oscillations corresponding to Fig.[71 other parameters 
are the same as in ^™Yb— ^^^Yb case. The panels and the lines arc the same as in Fig. [5l 
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FIG. 9: Time evolution of the quadrupole oscillations in the ^™Yb— ^^'^Yb mixture with the in-phase 
initial condition A_b=Af = 0.1. The upper and lower panels are Qb.f for the bosons and fermions, and 
the solid and dashed lines represent the results of TDGP + Vlasov and RPA calculations. 
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FIG. 10: Time evolution of the quadrupole oscillations in the ^^''Yb— ^^^Yb mixture with the out-of-phase 
initial condition \b = —^f — 0.1. The panels and the lines are the same as in Fig. [9l 




FIG. 11: Time evolution of Rz (solid lines) and Rt (dotted lines) for the quadrupole oscillation in 
i70y]3_i73y|3 mixture with the in-phase initial condition (corresponding to Fig. [9]); the upper and lower 
panels are for the bosons and fermions. 
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FIG. 12: The strength functions of the quadrupole oscillations with the in-phasc initial condition in 
170y13_i73y13 mixture (corresponding to Fig.[9|); the panels and lines are the same as in Fig. [5l 
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FIG. 13: The strength functions of the quadrupole oscillations with the out-of-phase initial condition in 
i70y|3_i73y]3 niixture (corresponding to Fig. [T0|: The panels and the lines are the same with Fig. [5l 
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FIG. 14: The fermion absolute strength functions of the quadrupole osciUations in ^^°Yb— "'^^^Yb gas with 
the in-phase and the out-of-phase initial conditions. 




FIG. 15: The absolute strength functions of the quadrupole oscillations in ^^"Yb— ^^^Yb mixtures with 
the in-phase initial conditions; the left (B1,B2,B3) and right (F1,F2,F3) panels are for the boson and 
fermion oscillations. The strength functions are deduced from the evolution in < t < 25 (the top 
panels), in 25 < t < 50 (the middle panels) and in 50 < r < 100 (the bottom panels), respectively. 
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FIG. 16: The absolute strength functions of the monopole oscillations in '^^''Yb— ^^^Yb mixtures with 
the in-phase initial conditions. The panels are the same as in Fig. 1151 




FIG. 17: Time evolution of the quadrupole oscillations in the ^^^Yb— ^^^Yb mixture with the in-phase 
initial condition = Xp = 0.1. The upper and lower panels are Qb.f for the bosons and fermions, and 
the solid and dashed lines represent the results of TDGP + Vlasov and RPA calculations. 
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FIG. 18: Time evolution of the qiiadrupole oscillations in the ^^^Yb— ^^^Yb mixture with the out-of-phase 
initial condition = — = 0.1. The panels and the lines are the same as in Fig. [9l 
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FIG. 19: Time evolution of the quadrupole oscillations in the ^^^Yb— ^^^Yb mixture with the in-phase 
(a) and the out-of-phase (b) initial conditions in the time-interval 60 < t < 120; the dashed and solid 
lines are for the boson and fermion oscillations, respectively. 



35 



4 - 
3 

« 2 
1 



-JwJ^A'^ 



CO ■Arty^'wv^- 



-1 - 



n — ^ — I — ' — r 
(a) 

TDGP + Vlasov 



■iwifWMPf|f||!|| 



(C) 
RPA 



||||jlj1l^Vljjn.3riAnj1JVliVI.«iArnrvWh^^ 




™iwniuyi|| 




(d) 



-10 



1.0 1.2 1.4 1.6 1.8 2.0 2.2 2.4.0 1.2 1.4 1.6 1.8 2.0 2.2 2.4 

FIG. 20: The strength functions of the quadrupole oscillations with the in-phasc initial condition in 
-'^^''Yb— ^^■^Yb mixture (corresponding to Fig. [T7)) : the panels and lines are the same as in Fig. [SJ 
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FIG. 21: The strength functions of the quadrupole oscillations with the out-of-phase initial condition in 
170y}5_i73y13 mixture (corresponding to Fig. [T8|: The panels and the lines are the same with Fig. \E\ 
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FIG. 22: The strength functions of the quadrupole osciUations with the in-phasc initial condition in BF 
mixture with qbe — 5.5976 (other parameters are the same as ^^''Yb— ^^'^Yb mixture); the panels and 
lines are the same as in Fig. [5] 
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FIG. 23: Time evolution of the quadrupole fermion oscillations in I'^^Yb-i^iYb (a), i^Yb-i'^^Yb (b) 
and ^^"Yb— ^^^Yb (c), where the boson motions are frozen. The initial condition is given by A_f — 0.1 for 
all cases. 
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FIG. 24: The fcrmion strength functions of the fermion quadrupole osciUations correspondmg to the cases 
in Fig. El 
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FIG. 25: Intrinsic frequencies of the quadrupole oscillations in boson (diamonds) and fermion (circles) 
components of the BF mixtures in the TDGP+Vlasov (full diamonds and circles) and RPA (open dia- 
monds and circles) calculations. The solid and dashed lines are for the sum-rule results of the boson and 
fermion quadrupole modes. 



